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Abstract Neonates with hypoxic-ischaemic (HI) brain injury were monitored 
using a broadband near-infrared spectroscopy (NIRS) system in the neonatal inten- 
sive care unit. The aim of this work is to use the NIRS cerebral oxygenation data 
(HbD = oxygenated-haemoglobin — deoxygenated-haemoglobin) combined with 
arterial saturation (SaO,) from pulse oximetry to calculate cerebral blood flow 
(CBF) based on the oxygen swing method, during spontaneous desaturation epi- 
sodes. The method is based on Fick’s principle and uses HbD as a tracer; when a 
sudden change in SaO, occurs, the change in HbD represents a change in tracer 
concentration, and thus it is possible to estimate CBF. CBF was successfully calcu- 
lated with broadband NIRS in 11 HIE infants (3 with severe injury) for 70 oxygen- 
ation events on the day of birth. The average CBF was 18.0 + 12.7 ml 100 g~! min"! 
with a range of 4 ml 100 g~! min“! to 60 ml 100 g~! min“. For infants with severe 
HIE (as determined by magnetic resonance spectroscopy) CBF was significantly 
lower (p = 0.038, d = 1.35) than those with moderate HIE on the day of birth. 
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1 Introduction 


Hypoxic ischaemic encephalopathy (HIE), which affects 1—2 full term newborns 
per 1000, is associated with severe neurodevelopmental problems and is the third 
leading cause of global child deaths [1]. HIE is an evolving process; after ischaemia 
during the injury, the brain is reperfused but the cerebral blood flow (CBF) can 
fluctuate between hypo-and hyper-perfusion in the first days after birth, and in some 
cases the injury can progress into a secondary energy failure [1]. Strategies to treat 
HIE focus on the maintenance of cerebral perfusion and preventing the cascade of 
events leading to further brain injury. 

We have previously presented a broadband NIRS device called CYRIL 
(CYtochrome Research Instrument and appLication) to monitor cerebral haemody- 
namics and metabolism via cytochrome-c-oxidase in the newborn brain [2]. In this 
paper, we present a method to additionally monitor CBF changes with CYRIL in 
neonates with HIE. Previously, researchers have used NIRS to measure CBF by 
using the oxygen swing method; oxygenation changes in the brain can be used to 
calculate the CBF via Fick’s law [3]. This states that the quantity of substance taken 
up by an organ per unit time (dQ/dt) is equal to the rate of arrival of the substance 
to an organ (arterial concentration (C,) multiplied by the blood flow (F)), minus the 
rate of departure from the organ (venous concentration (C,) multiplied by F): 


dQ/dt=F(C,-C,) (1) 


If Q varies with time, then this equation can be integrated with respect to time to 
obtain blood flow in terms of the quantity of tracer within the organ, and its arterial 
and venous concentrations: 


F=@Q/](C,-C,)dt (2) 
Q is equivalent to the product of the concentration of tracer in the tissue (Cr) and the 
unit weight of the tissue (W). This relationship can be used to obtain CBF in the 
more commonly used units, ml per 100 g of brain tissue per minute. Furthermore, if 
the accumulation is measured within the minimum transit time of the organ, reported 
as around 8 s for the brain in neonates [4], then the venous concentration of the 
substance will be zero. CBF can then be calculated from the ratio of the tracer accu- 
mulated to the quantity of tracer introduced during time t: 


F/W =C, /\C, (t)dt = CBF (3) 
When applied to NIRS measurements of CBF, the arterial concentration of sub- 
stance introduced over time t can be represented by the integral of the change in 
fractional arterial oxygen saturation (SaO,) multiplied by the total concentration of 
haemoglobin in the blood. The accumulation of substance in the organ becomes the 
concentration of oxy-haemoglobin (HbO,) in the brain, measured by NIRS, and 
CBF is then obtained from the equation: 
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CBF = K - A| HbO, | /| tHb-10° |f ASa0,dt 
with 
K = MWHb-10° / Dt -10 (5) 


where MW,, is the molecular weight of haemoglobin (64,500 g mol-!), tHb is 
the total haemoglobin concentration, and Dt is the tissue density, assumed to be 
1.05 g.ml~! in the brain [3]. 

Haemoglobin difference (HbD = HbO, - HHb) is a measure of oxygen delivery 
that has a larger amplitude (and thus higher signal to noise) and is also more stable 
than either of the two individual signals as artefacts are subtracted. If total haemo- 
globin concentration remains stable throughout the measurement, the change in 
HbD can be expressed as twice the change in HbO, and therefore CBF is: 


CBF = K-A[HbD]/2| tHb-10° |f ASaO, dt is 

Here we demonstrate the use of this method not with controlled oxygen chal- 
lenges, but on spontaneous desaturation events in newborns with brain injury under- 
going intensive care. This allows non-invasive measurements of CBF at the cot side 
without intervention. 


2 Methods 


2.1 Data Collection 


Ethical approval for the Baby Brain Study at University College London Hospitals 
(UCLH) Trust, London was obtained from the North-West Research Ethics Centre 
(REC reference: 13/LO/0106). Term infants born at or transferred to UCLH for 
treatment of acute brain injury were eligible for investigation; infants with congeni- 
tal malformations were excluded. Informed consent was obtained before each study. 
Each subject was monitored continuously with EEG and treated with hypothermia 
which was initiated within 6 hours of birth; body temperature was lowered to 
33.5 °C and was maintained for 72 hours before gradual rewarming. 

Broadband NIRS measurements were collected with CYRIL continuously over a 
period ranging from 1—15 hours on the first day of life during therapeutic hypothermia. 
NIRS probes were placed over the left-side of the forehead (covering the left frontal 
lobe) with 3 cm optode separation and with a sampling frequency of 1 Hz. Changes in 
chromophore concentrations (HbO,, HHb, and the oxidation state of cytochrome-c- 
oxidase, oxCCO) were calculated from the measured changes in light attenuation using 
the modified Beer-Lambert law. This was applied across 770-906 nm with a differential 
pathlength factor of 4.99 [5], which is corrected for the pathlength wavelength-depen- 
dency. Examples of obtained intensity spectra are published in Bale et al. [2]. 
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Systemic data from the Intellivue Monitors (Philips Healthcare, UK) were col- 
lected using ixTrend (ixellence GmbH, Germany). Signals recorded include oxygen 
saturation (SaO,) measured by pulse oximetry, heart rate (HR), respiratory rate 
(RR), and mean arterial blood pressure (MABP). 

An MRI scan was performed after rewarming; this occurred on average on day 7 
of life (range: day 5—15). MR scans included a measurement of thalamic Lac/NAA 
with proton ('H) magnetic resonance spectroscopy (MRS); Lac/NAA score > 0.3 
has been shown to be associated with poor outcome [6]. 


2.2 Data Analysis 


Data analysis was carried out in MATLAB (Mathworks, USA). Systemic data were 
down-sampled and interpolated to the NIRS data timeframe (1 Hz). Changes in 
SaO, of 2% or greater occurring in under 7 seconds (determined to be the transit 
time of the brain) were identified manually using a graphical user interface written 
in-house. Events were rejected if the initial SaO, was less than 70%, if there was a 
change in HR or RR of more than 10 bpm during the event, or if the HbD standard 
deviation was greater than 1 uM in the 10s before the event. An example of the pulse 
oximetry and broadband NIRS data used is shown in Fig. 1a and b). The gradient of 
the change in HbD and integral of SaO, with respect to time during the first 7 sec- 
onds of the event was used to calculate the CBF (see Fig. 1c) and Eq. 6). Data are 
presented as median and confidence interval, and statistical significance was 
assessed with the Mann-Whitney U test (p > 0.05 is significant) and the effect size 
was tested by Cohen’s d (d > 1.2 is ‘very large’) [7]. 
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Fig. 1 Example of the CBF calculation method: (a) SaO, and (b) HbD. The vertical lines indicate 
the region of interest (10 seconds) used for the calculation. (c) The change in HbD plotted against 
the equivalent cumulative change in SaO, integral; with a fitted fourth order polynomial and a fit- 
ted linear regression. The CBF in this event is 22.1 ml 100 g~' min 
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3 Results 


CBF was calculated during 42 desaturation episodes noted in 11 infants with HIE on 
the day of birth. Of these infants, 3 were predicted to have a poor outcome by the 
MRS Lac/NAA biomarker. See Table 1 for clinical details, number of events and 
mean CBF for each neonate. The mean gestational age was 39.3 + 1.2 weeks and the 
mean birth weight was 3262 + 334 g. The median CBF was 12 ml 100 g~! min“! with 
a range of 4 ml 100 g~! min"! to 60 ml 100 g"! min™!. For infants with severe HIE, 
CBF was lower (8.5 (6.5—11.0) ml 100 g~! min~!) than moderate HIE infants (15.2 
(4.0-60.0) ml 100 g~! min=!) on the day of birth (p = 0.038, d = 1.35, see Fig. 2). 


4 Discussion 


The NIRS measured CBFs are in line with those previously measured in newborns 
by the oxygen swing technique [4, 8], Xenon clearance [4, 9], dynamic contrast- 
enhanced NIRS with indocynanine green [10] and phase contrast MRI [11] which 
had means of 29.7 (range 11.4-62.1), 11.7 (range 5—26.1), 14.8 (range 5.1—30.4), 
13 + 7.4, 16.5 + 2.1 and 9.6 (range 7.5—15.1) ml 100 g~' min“! respectively. This 
suggests that using the oxygen swing method on spontaneous desaturation events is 
a feasible way to measure CBF without intervention. 

Neonates with severe HIE experienced hypoperfusion on the first day of life. This 
is consistent with a finding from an MRI arterial spin labelling study; infants with 
severe HIE had significantly lower CBF on the first two days of life than those with 
moderate HIE [12]. The decreased CBF in more severe injury could be caused by 
evolving pathophysiology of HIE - following HI, there is a period of hypoperfusion 


Table 1 Clinical details and number of eligible desaturation events for neonates included. 
Unfavourable outcome was determined by MRS Lac/NAA > 0.3. Birth weight measured in grams, 
CBF measured in ml 100 g™! min“! 


GA Birth 
(weeks) 


| Injury 
outcome 


Neonate | Gender Lac/NAA No. of events 


F 0.2 Favourable 2 14+4 
M 0.16 Favourable 6 10+1 
F 0.08 Favourable 12 16+5 
M 0.25 Favourable 4 47+ 10 
F 0.15 Favourable 9 11+4 
F 0.35 Unfavourable | 2 71 

F 0.43 Unfavourable |4 9+2 

F 0.14 Favourable 4 21411 
F 0.2 Favourable 2 41+2 
M 0.34 Unfavourable |1 11+0 
F 0.15 Favourable 3 12+3 
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Fig. 2, Boxplot of CBF on first day of life per injury severity group (*p = 0.024) 


associated with hypometabolism during the latent phase, followed by relative hyper- 
perfusion in the secondary phase [1]. 

Therapeutic hypothermia itself can reduce CBF; an index of CBF measured by 
diffuse correlation spectroscopy has been shown to be reduced during cooling com- 
pared to after rewarming in a HIE group, as well as compared to a control group at 
normothermia [13]. However, this cannot account for the differences in CBF 
between the injury severity groups because all neonates were being treated with 
hypothermia during the measurement. It important to be aware of the impact of 
cooling on CBF, as an early study without hypothermia showed the opposite trend; 
there was a trend for higher CBF in infants with severe HIE [8]. This observation of 
hypoperfusion solely in the severe group suggests that the temporal progression of 
injury and/or the response to treatment is injury dependent. 

We have demonstrated a method of calculating absolute CBF non-invasively and 
without any intervention using NIRS during spontaneous hypoxic events in HIE 
infants. This technique is advantageous because it can be performed at the cot side 
and does not require any intervention, but is limited because the spontaneous nature 
of the events. Indeed, desaturation episodes are not guaranteed in all infants during 
a study period. However, we could identify a difference in the level of perfusion 
between neonates with moderate and severe injury on the first day of life. This sug- 
gests that monitoring CBF during HIE will give important clinical information that 
may have further prognostic value. 
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